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Abstract: 

The transition from the near to the far field of the fundamental mode 
radiating out of a photonic crystal fiber is investigated experimentally 
and theoretically. It is observed that the hexagonal shape of the near 
field rotates two times by 7r/6 when moving into the far field, and 
eventually six satellites form around a nearly gaussian far-field pattern. 
A semi-empirical model is proposed, based on describing the near field 
as a sum of seven gaussian distributions, which qualitatively explains 
all the observed phenomena and quantitatively predicts the relative 
intensity of the six satellites in the far field. 
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1 Introduction 

Photonic crystal fibers (PCF) are a new class of optical fibers which has revealed many 
surprising phenomena and also holds a big promise for future applications (see e.g. [1, 2, 




Fig. 1. Schematic of a single-mode PCF (z < 0) with an end-facet from where light 
is radiated into free space (z > 0). 

3]). These PCFs are made from pure silica with a cladding consisting of a regular lattice 
of air-holes running along the fiber axis. Depending on the arrangement of the air-holes 
the guiding of light can be provided by either modified total internal reflection [4, 5] or 
by the photonic band-gap effect [6, 7] and PCFs can even be endlessly single-mode [8] 
because of the wavelength dependence of the cladding index. For the basic operation we 
refer to the review of Broeng et al. [9] . 

Understanding the shape and radiation pattern, as illustrated in Fig. 1, of the mode 
in the endlessly single-mode PCF is very important. E.g. in tests and applications this 
is essential for estimations of coupling efficiencies and for determining the mode field 
diameter from the far-field distribution. Furthermore, it is fundamentally the simplest 
structure with a hexagonal cladding, and hence the understanding of this structure will 
be a natural basis for understanding the modes of more sophisticated PCF structures. 
In this paper we present a semi-empirical model which is capable of explaining both 
the near and far-field distribution of the mode, but most importantly also accounts 
for the fine structure in the transition from the near to the far field. The simplicity of 
the model allows for a phenomcnological interpretation of the shapes of the near and 
far-field patterns. 

2 Experiments 

The measurements reported arc for a PCF with a triangular air-hole lattice with pitch 
of A ~ 3.5 /im and air holes of diameter d/A ~ 0.5. The measurements reported here 
were performed at a free-space wavelength of A = 635 nm, where the light is guided in 
a single mode in the silica core of the fiber formed by a "missing" air hole. In panel a 
of Fig. 2 a micro-graph of the fiber structure can be seen. 

The near-field distribution was measured using a microscope objective to magnify 
the mode onto a Si-based CCD camera. In Fig. 2b the intensity distribution is shown 
at focus. By translating the fiber away from the focal plane, the intensity distribution 
may be imaged at different distances between the near and the far field. This is shown 
in panels b to h in Fig. 2. As expected the mode at focus has a hexagonal shape, that 
extends in the six regions between the inner holes and is sharply confined at the six silica- 
hole interfaces. However, when the image is defocused, the shape at first transforms into 
a nearly circular shape (panel c) followed by a hexagonal shape rotated by an angle of 
7r/6 with respect to the focus (panel d). After this the shape again becomes close to 
circular (panel e) , and finally transforms into the original hexagonal orientation (panel 
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Fig. 2. Experimentally observed near-field intensity distributions for a PCF with 
A ~ 3.5/im and d/A ~ 0.5 (micro-graph in panel a) at a free-space wavelength 
A = 635 nm. The distance from the end-facet varies from z = to z ~ 10 (im 
(panels b to f). At a further distance the six low-intensity satellite spots develop 
(panels g and h, logarithmic scale). 

f) with six satellites emerging from the distribution (panels g and h). It is noted that 
the orientation of the satellites is rotated by ir/6 with respect to the six inner holes 
surrounding the core. In Fig. 3 (right) the intensity distribution in the far-field limit 
is shown (several centimeters from the fiber end- facet), obtained using a commercial 
far-field profiler. Here, the satellites have fully developed and as shown in the cross 
sectional plot in Fig. 3 (left) the peak intensities of the satellites are more than two 
orders of magnitude lower than the main peak. Hence, a reasonably accurate analysis 
of the far field may be performed considering only the main peak. 

Apart from being a fascinating and intriguing evolution of the mode shape from the 
near to the far field, it is important to be aware of these transitions in any application 
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Fig. 3. Experimentally observed far-field intensity distribution showing an overall 
gaussian profile with six additional low-intensity satellite spots along one of the two 
principal directions (line 2). Angles are given in radians. 
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Fig. 4. Panel a shows the experimentally observed near-field intensity along the two 
principal directions 1 and 2 (see insert of panel b). Panel b shows the numerically 
calculated intensity distribution in a corresponding ideal PCF with the solid lines 
showing the intensity along the principal directions and the difference. The blue 
and red dashed lines show gaussian fits to I2 and I2 — Ii and the dashed green line 
shows their difference. 



that involves imaging of the modes. E.g. for estimations of the mode field diameter 
and effective area based on near-field analysis, it is important to focus the mode cor- 
rectly, and the positions corresponding to panel b and panel d in Fig. 2 may easily be 
confused. They both show the hexagonal shape and have only slightly different mode 
sizes. Hence, as a measurement procedure for determining the mode field diameter, a 
direct measurement of the near field may be even more tricky than it is for "standard 
technology fibers" with circular symmetry. 

In panel a of Fig. 4 two cross-sections of the measured near-field distribution are 
shown, one taken along a line passing through opposite hole centers (1) and the second 
taken along a line passing between the holes (2) (rotated by an angle 7r/6 with respect 
to the first). It is noted that a gaussian distribution is a significantly better fit to the 
intensity along line (2), motivating a simple interpretation of the mode shape : the mode 
is a circular gaussian distribution from which a narrow distribution localized at each of 
the centers of the six inner holes is subtracted. This simple interpretation is theoretically 
modeled in the following. 

3 Theory 

In order to simulate the radiated field we start from the fully-vectorial fields in the 
single-mode all-dielectric PCF 



where h(x, y) and /3(u>) are the transverse fields and the propagation constant, respec- 
tively. These we calculate numerically by a fully-vectorial plane- wave method [10]. 

Substantial insight in the physics of the radiation problem can be gained by ex- 
panding h(x, y) in gaussians. Introducing the notation s = (x, y) and using that the 
components of h(x,y) can be chosen either real or imaginary we consider 



H(x,y,z) = h(x,y)e ± ^ z 
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For the radiation into free space this gives a linear combination of expanding gaussian 
beams and this is a well-studied problem, see e.g. [11, 12]. Neglecting the small back- 
scattering from the end-facet, the gaussian u(s,w) at finite z > transforms as 
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Fig. 5. Near-field intensity distribution calculated from Eq. (5) with values of u>h, 
w c , and 7 determined from the intensity in the PCF obtained by a fully- vectorial 
calculation, see Fig. 4. The distance varies from z = to z = 8A (panels a to 
i) in steps of Az = A (see also animation with Az = A/4, http: //www. crystal- 
fibre, com/technology /movie. gifi 3 Mbyte). 
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where R{z) = z(l + fc 2 w 4 /4z 2 ) and W(z) = w(l + 4z 2 /A; 2 w 4 ) 1 / 2 . In the following we 
consider a particular simple linear combination in the PCF; 



I(s) = A 2 u(s, w c ) — 7 u(s — sj 



Wh 



(4) 



where Sj = R c x (cos ^p, sin with (A/R c )sj being the center position of the six 
air holes nearest to the core. Here, R c ~ A (the radius of the silica core), w c ~ A (the 
mode-field radius), and Wh ~ d/2 (the radius of the air holes). The first term gives the 
over-all gaussian intensity profile of the mode and with 7 ~ u(R c , w c ) the additional six 
terms of opposite sign suppress the intensity at the six air-holes nearest to the core. For 
finite z > the intensity transforms as 




In panel b of Fig. 4 we show an example of the intensity distribution in an ideal PCF 
with A = 3.5 [im and d/A = 0.5 at A = 635 nm corresponding to experimental situation. 
For the dielectric function we have used e = 1 for the air holes and for the silica we 
have used e = (1.4572) 2 = 2.123 based on the Sellmeier formula. While Eq. (5) may 
seem too simplistic the good fits to gaussians strongly justify it and as we shall see it 
reproduces the physics observed experimentally. 

In Fig. 5 we show the corresponding near field based on Eq. (5). The profile at the 
end-facet (panel a) first transforms into a close-to-circular profile (panel c) followed by 
a hexagonal shape rotated by ir/6 (panels d to f), a close-to-circular profile (panel g), 
and finally a hexagonal shape (panels h and i) with the same orientation as at the end- 
facet (panel a) . Comparing with Fig. 2 this is qualitatively in excellent agreement with 
the experimental observations. The fact that the fully coherent scattering description 
qualitatively reproduces the experimentally observed tt/6 rotation gives strong indica- 
tions of its nature; it is a phenomena caused by an interference between the different 
gaussian elements used in the decomposition of the fundamental mode in the PCF. In 
Fig. 6 we show the corresponding intensity distribution in the far-field limit which is 
in a very good agreement with the experiments, see Fig. 3. It is seen that the satellites 
are reproduced and are in fact oriented in the same way as in the experiment. Moreover 
the relative intensities between the satellites and the main peak in Fig. 6 (left) are very 
similar to the ones in Fig. 3 (left). 

Finally, for fibers where the air holes modify the overall gaussian profile sufficiently 
(not shown) we find indication of additional higher-order spots further away from the 
center of the intensity distribution which can also be seen experimentally. 



4 Conclusion 

The evolution of the mode shape of a PCF with a triangular cladding has been inves- 
tigated in the transition from the near to the far field. When moving away from the 
near field at the focus of the fiber end- facet, it has been observed that the hexagonal 
orientation is rotated two times by 7r/6 after which six satellites emerge in the radia- 
tion pattern. In the far-field limit the satellites remain in the pattern, having a relative 
peak intensity more than two orders of magnitude less than the main peak and with 
an orientation rotated by 7r/6 relative to the six inner holes around the fiber core. All 
these observations have been reproduced theoretically, by approximating the near-field 
distribution by a main gaussian peak from which six narrow gaussians located near the 
center of the six inner holes have been subtracted. From the simulations it is concluded 
that the changes of shape in the radiation pattern are caused by an interference between 
the different gaussian elements used in the decomposition of the fundamental mode in 
the PCF. 

The results presented here are very important for understanding and analysing the 
behavior of the mode in many optical systems based on photonic crystal fibers - es- 
pecially those involving imaging and focusing the mode. Furthermore, the successful 
idea of decomposing the near field of the mode into seven localized distributions can 
be adapted in future work aimed at a simple quantitative description of the near and 
far-field distributions, particularly for relating the measured far field to the physical 
structure around the fiber core that influences the near field. The latter is very in- 
teresting in the field of fiber measurement procedures, since the far-field analysis of 
conventional optical fibers cannot directly be adapted to PCFs because of the lack of 
cylindrical symmetry 
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